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<L> 

■ The scattering of light at large wave-vector magnitudes k in nematic systems containing topo- 

logical defects is investigated theoretically. At large k the structure factor S(k) is dominated by 
power-law contributions originating from singular order-parameter variations associated with topo- 
logical defects and from transverse thermal fluctuations of the nematic director. These defects 
(nematic disclinations and hedgehogs) lead to contributions of the form pAk~^ ("the Porod tail"), 
where p is the number density of a given type of defect, A is a dimensionless Porod amplitude, and 
£ is an integer- valued Porod exponent. The Porod amplitudes and exponents are calculated for all 
types of topologically stable defects occurring in uniaxial and biaxial nematics in two or three spatial 
dimensions. The range of wave-vectors in which the contributions to the scattering intensity due to 
defects dominate the contribution due to thermal fluctuations is estimated, and it is concluded that 
for experimentally accessible defect densities the range of observability of the Porod tail extends 
over one to three decades in scattering wave- vector magnitude k. Available experimental results on 
phase ordering in uniaxial nematics are analyzed, and applications of our results are suggested for 
q ■ light-scattering studies of other nematic systems containing numerous defects. 
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PACS numbers: 61.30.Jf,78.20.-e,64.60Cn,82.20Mj 
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I. INTRODUCTION 
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CN| ■ When a liquid crystalline material is brought from the isotropic phase to the nematic phase, the turbidity (i.e., the 
total intensity of scattered light) of the material increases by a factor of the order of 10 6 , and the sample appears 
cloudy [jlj. In samples in which the nematic ordering is well aligned, this intense scattering is due predominantly to 
thermal fluctuations in the orientation of the nematic director, and has been extensively studied and reported in the 
literature jjlp],plB] . If, on the other hand, the sample is not well aligned and contains a large number of topological 
defects (disclinations, nematic hedgehogs, or surface defects), the turbidity becomes dominated by scattering from 
essentially static director inhomogeneities associated with these defects pj. 

A particularly visible manifestation of defect-associated scattering is the so-called "Porod tail" contribution to the 
scattering intensity. For sufficiently large magnitudes k of the scattering wave-vector k, the scattered light intensity 
O ■ /(k) of a nematic system containing topological defects decays as a power law k~~ x , where \ is an integer- valued 
exponent that depends on the dominant type of defect present in the system. Such behavior was confirmed and the 
values of \ have been discussed in both the experimental ||U and the theoretical literature on the kinetics of phase 
ordering of nematic systems containing numerous defects generated during a quench from the isotropic to the nematic 
phase. A full calculation of the form of the power-law tails associated with nematic defects, and an evaluation of 
■ the relative importance of these contributions compared to the fc~ 2 contribution associated with thermal fluctuations 
have, however, not yet been given in the literature, and form the subject matter of the present Paper. 

Our results are, in principle, applicable to any nematic system containing topological defects, but should be es- 
pecially useful in interpreting detailed measurements of scattered-light intensity in highly disordered systems (i.e., 
those containing numerous defects). Such configurations arise in low- molecular- weight nematics for example when the 
transition from the isotropic to the nematic phase is sudden (e.g., induced by a temperature quench) |p|,p|,p|J^, ^0| , pd[ | , 
or when an ori gina lly well-aligned nematic sample is put in a sufficiently strong shear flow JT2fl or under high alter- 
nating voltage |L3|]. In polymer nematics, numerous defects often persist even in the absence of any external agents 
|i~4].[l5|], and significantly affect the mechanical jl6| and electro-optical |l7j] properties of these materials. Our results 
permit the extraction of information on the type and number of topological defects present in such systems and the 
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monitoring of the dynamics of the defects. As our results are exact in the appropriate scattering wave-vector range, 
they can also be used to test the validity of analytical theories and simulations of phase ordering in nematics. 

This Paper is organized as follows. In Sec. || we discuss the general issue of the origins of power-law contributions 
to the structure factor S(k) (i.e. the Fourier-transformed order-parameter correlation function) in systems possessing 
continuous symmetries. In general, two types of power-law contributions of quite distinct origins are present: those 
due to transverse thermal fluctuations of the order parameter, and those due to topological defects. The order- 
parameter variations associated with topological defects give rise, at sufficiently large k, to contributions of the form 
pAk~^ (the Porod tail), where p is the number density of a given type of defect, A is a dimensionless Porod amplitude 



and £ is an integer- valued Porod exponent. In Sec. Ill we re-derive and generalize some of the results of Bray and 
Humayun |l8f| for Porod amplitudes and exponents in the O(N) vector model by using a somewhat different and 
computationally simpler method. In the remainder of the Paper, we use this method to calculate Porod amplitudes 
and exponents for topological defects in uniaxial and biaxial nematic liquid crystals. The forms of the Porod tail that 
correspond to hedgehog defects, disclination lines, and ring defects in uniaxial nematics are calculated in subsections 



fVB| - |lV4 A case of special interest is that of a wedge-type ring defect (i.e. disclination loop) in a uniaxial nematic, 



where two separate Porod regimes (having distinct exponents and amplitudes) arise for length scales larger than and 
smaller than the radius of the disclination loop. The presence of two Porod regime s is s pecific to the nematic case, 
and does not occur for any type of defect in O(N) vector model systems. In Sec. IV Fj we discuss the influen ce of 



transverse thermal fluctuat ions of the nematic director on the large-A: behavior of the structure factor. In Sec. IV F 
we use the results of Sees. I V A - 1 V E to analyze the results of the light-scattering experiments of Refs. H[|, which 



investigate the process of phase-ordering kinetics in uniaxial nematics. In Sec. ^ we generalize our results for the 
uniaxial nematic to the case of non-abelian defects in biaxial nematics, the dynamics of which have recently been 
investigated experimentally [|l9| and theoretically p0| . A general discussion of corrections to the results derived in 
this Paper that may arise from effects such as defect interactions, defect curvature, and the presence of the defect 



analysis of experimental data. 



core, is given in Sec. VI. We conclude, in Sec. VII, with a summary of our results and suggestions for their use in the 



II. GENERALIZED POROD LAW 



Consider an ordered system characterized by an order-parameter field $(r), where 3? is an A-component vector in 
order parameter space and r is the radius-vector in d-dimensional real space. The principal quantity of interest in 
many theoretical and experimental investigations of ordered systems is the structure factor 5(k), which is the Fourier 
transform of the real-space correlation function C(r) of the order parameter. This quantity, which is directly measur- 
able via the appropriate scattering experiment, characterizes the degree of inhomogeneity of the order parameter at 
length scales of order fc _1 . By definition, 

5(k) ee J d d r e lkr C(r) , (1) 
with the real-space correlation function C(r) given by 

C(r) = -J— /d^*(x).*(x + r), (2) 

M 0(N) J 

where the integration is taken over the whole system. Here A/q(n) is a normalization factor, 

M 0( N) = J d d x*(x)-*(x), (3) 

chosen to ensure that C(r)| r= o = 1. Notice that with this choice of normalization C(r) is dimensionless, whereas S(k) 
has the dimensions of a volume. Equivalently the structure factor may be expressed as 

5(k) = TT^*(k) • *(-k) , (4) 

Mo{N) 

where $(k) is the Fourier-transformed order parameter, 

*(k) ee J d d re ik r $(r). (5) 
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For systems that are isotropic on macroscopic scales (such as a bulk system undergoing phase ordering), 5(k) depends 
on k through the magnitude k = |k| only. 

We now discuss the contributions to S(k) that have the form of a power law ( i.e., ak~ x ) where a and \ are 
k- independent constants. First, we briefly recall the effect of thermal fluctuations. Consider an O(N) vector model 
system (with N > 2) in the ordered phase with the thermally-averaged value of the order parameter <fr(x) pointing in 
the x\ direction. Due to the O(N) symmetry, fluctuations in ^(x) perpendicular to the x\ direction (i.e., transverse 
fluctuations) cost an arbitrarily low amount of free energy for any large- length-scale fluctuations; this leads to strong 
scattering at small k in the whole temperature range of the ordered phase. [In contrast, strong low-fc scattering in 
systems having a scalar order parameter (for which transverse fluctuations do not exist) arises only in the immediate 
vicinity of a critical point ("critical opalescence") and is due to longitudinal fluctuations (i.e., fluctuations in the order 
parameter magnitude).] For our purposes, the important property of transverse fluctuations is that they occur on all 
length scales between the size of the system (small k) and the microscopic coherence length of the order parameter 
(large k). In a system of spatial dimension d > 2 |22j and described by the standard form of the gradient free energy E 



K 



E=~J ^arV^xjV^tx) (6) 

(where summation over the indices i and j is implied), the contribution to the structure factor due to thermal 
fluctuations has the form a(T)k~ 2 with a temperature-dependent amplitude given by the microscopic length scale 
a(T) = k B T/K @. 

Next, we discuss power-law contributions to the structure factor that arise from topological defects. Topologically 
stable defects occur in all O(N) vector model systems whenever d > N. The order-parameter variations associated 
with the defects give rise to the Porod-tail form of the structure factor, S(k) ~ k~ x , for sufficiently large values of k. 
In contrast to the thermal-fluctuation contribution, the power-law contributions to S(k) due to topological defects do 
not depend on temperature (provided the defect structure is not affected by the temperature) and occur in systems 
with discrete symmetries as well as in those with continuous symmetries. 

In scalar systems, this power-law behavior has long been understood as arising from abrupt changes in the order 
parameter across domain walls present in the system |p4[p5| . Consider, for simplicity, a one-dimensional system in 
which the order parameter suffers a kink. Viewed on length scales larger than the kink width (which is given by the 
coherence length of the order parameter), the order parameter abruptly changes from —I to +1 at the kink location. 
Consequently, the integrand in Eq. (|^) equals — I within the distance \r\ from the kink (domain wall) location and +1 
elsewhere in the system. This results in the normalized real-space correlation function of the form C(r) = 1 — 2\r\. 
Notice that C(r) is non- analytic at r = 0, reflecting the singular variation of the order parameter at the kink. The 
non-analytic part of C(r), after Fourier-transforming, yields a power-law form of the structure factor, S(k) ~ A; -1 , 
valid for values of k larger than the inverse system size and smaller than the inverse kink width. 

For domain walls in dimensions d higher then one, additional factors of 1/fc arise due to the spatial extent of the 
domain walls, with the result (averaged over domain- wall orientations) S(k) cx k~^ d+1 \ In regions between the domain 
walls, the order parameter does not vary, and therefore there are no extra contributions to the structure factor for 
k values larger than the inverse domain size. We thus obtain a structure factor that has a power-law form for k in 
the range between the inverse average domain wall separation and the inverse domain wall width, with the pre-factor 
proportional to the total domain wall area in the system. 

Bray and Humayun |fl8f investigated in detail how the singular variation of the order parameter around a point or 
line defect in a system with continuous symmetry will likewise lead to an asymptotic power-law dependence on k of 
the structure factor S{k). In Ref. |p^| , they calculated exactly the exponent and the amplitude in the contribution 
Ak~ x to the structure factor coming from an isolated defect order-parameter configuration in the O(N) vector model. 
They also argued that for sufficiently large k, the structure factor of a system containing a number of topological 
defects can be obtained by multiplying Ak~ x by the number density of defects. The reason is as follows. Variations in 
the order parameter on very short scales [probed by S(k) with k large] occur only in the singular regions surrounding 
the defect cores, which are separated by regions that barely contribute to S(k). It is therefore appropriate to treat the 
contribution from each defect independently, and as coming from an isolated defect, provided that the configuration 
in the vicinity of the defect core is not affected by the presence of other defects. While this assumption is known to 
be valid in some special cases [such as the 0{2) vector model], it is in general non-trivial. For example, it was argued 
in Ref. j2(| that, for the case of point defects in the 0(3) vector model in d = 3, a highly asymmetrical, string- like 
configuration develops in between the defects. Nevertheless, the prediction for S(k) based on the assumption of 
undeformed defect configurations has been found to be in good agreement with num erica l simulations in several O(N) 
vector model phase-ordering systems 27 . In the calculations we present in Sees. III-V|, we shall assume that the 



regions close to the defect cores in a phase ordering system do possess the structure of an isolated defect, and in 



Sec. VI we shall comment on the nature of possible corrections to our results. 
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The discussion in the preceding paragraphs leads us to the following formulation of the generalized Porod law for 
the orientation-averaged structure factor at zero temperature of a system containing topological defects: 



S(k) ~ Apk~ x , (7) 

where p is the defect density (i.e., the domain wall area, string length, or number of point defects, per unit volume of 
the system), and A is a dimensionless amplitude characterizing the given type of defect. The generalized Porod law 
is expected to be valid in the range L _1 <C k -C £ _1 , where L is the average separation between defects present in the 
system and £ is the size of the defect cores. Specific cases exemplifying the validity of Eq. ((?]) will be given throughout 
the rest of the paper. The value of the Porod exponent \ can be anticipated on purely dimensional grounds. Denoting 
by d the spatial dimensionality of the system and by D the dimensionality of the topological defect (that is, D = 
for point defects, D = 1 for string defects, and D = 2 for domain walls), we see that the defect density p in Eq. (0) 
scales as L D /L d , whereas the structure factor S(k) has the dimension of L d . Consequently, the exponent \ is given 



simply by 2d — D (see Sec. VI for a more precise derivation of this result). In contrast, the dimensionless amplitude 
A depends in more detail on the order parameter in question. For example, we shall see that the amplitudes A for 
the nematic hedgehog and for the 0(3) monopole in d = 3 differ substantially. 

Nematic liquid crystals represent an interesting case where several distinct power-law contributions to the str uctur e 



factor can be present simultaneously in the same system. These contributions are discussed in detail in Sees. IV B- 



IV E and are due to line defects (disclinations), point defects (nematic hedgehogs) and thermal fluctuations. The 
competition among these terms can lead to a variety of forms of the tail of the structure factor appropriate to 
distinct physical conditions. Before describing in detail the case of nematic liquid crystals, however, we illustrate our 
computational techniques on the case of Porod amplitudes and exponents in 0(V)-symmetric vector-model systems. 



III. O(N) SYMMETRIC VECTOR SYSTEMS 

The Porod exponent x arL d amplitude A of a contribution to <S(k) from a single defect in the configuration of an 
O(N) vector model was computed in Ref. by first calculating the real-space correlation function C(r), Eq. (||), and 
then Fourier-transforming the result to obtain S(k). During this calculation, certain analytic terms (see Ref. ]l8|), 
which do not contribute to the power-law part of the Fourier transform, were discarded; also, for O(N) systems 
with N even, the calculation was complicated by the appearance of logarithmic terms in the intermediate result for 
C(r). We shall show below that the results of Ref. [Q can be obtained in a more simple way by working directly in 
Fourier space, i.e., by first Fourier-transforming the order-parameter configuration, and then using Eq. (|j) to obtain 
the structure factor |p8| . Besides being shorter, this method does not discard any terms, and does not involve any 
logarithmic terms. 



A. Point defects in O(N) symmetric vector systems 

We first perform the calculation for the case of the O(N) vector model system in d = N spatial dimensions. In 
this case, only point defects (i.e., monopoles) are topologically stable; moreover, for d — N = 2, monopoles with 
charge amplitude \Q\ > 1 are energetically unstable towards the decay into monopoles with charge +1 or —1. In the 
following, we consider the case Q = +1; an identical result for the structure factor is obtained in the case Q = — 1. 

The configuration of the Q = 1 (radial) monopole in a system with the free energy given by Eq. (^J) is described by 

*(r) = - . (8) 

r 

Note that we have omitted any spatial variation in the amplitude |4>| of the order parameter. This corresponds 
to our focusing on length scales larger that the so-called core of the defect (i.e., the region of space in which the 
competition between condensation and gradient free energies gives rise to significant modification of |3?|). For a 
configuration having constant the correlation function Eq. (Q) and the structure factor Eq. (^) do not depend on 
|«&|; consequently, we have chosen |<&| = 1 for our calculation. 
Due to spherical symmetry, the Fourier transform of <&(r) takes the form 

*(k)=k/(fc 2 ), (9) 

where / is a certain scalar function of the scalar k 2 . By taking the scalar product of Eq. (^|) with k and solving for / 
we see that we may write 
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S(k) = k k- 2 J d d r r- 1 k • r e ak r = -ik fc" 



F{\ 2 k 2 ), (10) 

A=l 



where the function F is defined via 

F{\ 2 k 2 ) = J d d rr- x e iX * r , (11) 

i.e., F is the d-dimensional Fourier transform of the Coulomb potential r _1 , which can readily be shown |2l)| , for 
k ^ 0, to be 

( 47r )(d-i)/2 / d-1 



W^) = ^3=i-H— J' (12) 

where T(x) is the standard Gamma-function. By evaluating the derivative of F, and inserting it into Eq. (|l0|), we 
find 

*(k) = i(d - i)(^r~^ 2 r A,. (is) 

The normalization factor Mq(n)j given by Eq. (||), is the volume V of the system. It follows from Eq. (||) that the 
structure factor S(k) of a single point defect is then given by 

S(*)-ii(^f(i±l)^. (14) 

Specifically, we obtain S"(fc) = y _1 47r 2 fc -4 for an 0(2) vortex in d = 2, and = V~ 1 Ylic'k~^ for an (9(3) monopole 
in d = 3. Note that the result Eq. ( |T4| ) does not depend on the magnitude 3>| of the order parameter (we have chosen 
1 4? | = 1 in this subsection) as both the normalization factor Eq. (g|) and the numerator in Eq. (|]) are proportional to 

For a system having p defects per unit volume, the normalization factor factor V 1 is replaced by p. The expression 
Eq. (jlj) is then in agreement with the result for S(k) given in Ref. fl8|] . It should be noted that by fully exploiting the 
spherical symmetry at hand, the method used above enables us to pass from a calculation involving vector quantities 
to a calculation involving only scalars. In Sec. IV B| , we shall similarly be able to pass from a tensorial calculation, 
appropriate for the case of the nematic order parameter, to a scalar one. 



B. Vortex lines in three-dimensional 0(2) vector systems 

We now turn to the case of line defects in O(N) vector systems, specifically concentrating on the physically prominent 
case of vortex lines in three-dimensional 0(2) vector systems. 

In the case of the spherically-symmetric point defect in the O(N) vector system, investigated in the previous 
subsection, the final result, Eq. (|14|), is independent of the orientation of the scattering wave- vector k, due to the 
O(N) symmetry. In contrast, the structure factor of a segment of a line defect depends on the angle between k and 
the orientation of the line segment; specifically, S(k) is dominated by contributions from defect segments that are 
perpendicular to the vector k, as there are no short-distance order-parameter variations in the direction along the 
defect line. 

We first discuss the structure factor of a single straight vortex line. Consider a vortex-line segment in d = 3 having 
its core located on the line (x, y) — (0,0), and extending from z = —L/2 to z = L/2. The order parameter $ does 
not depend on z, and the minimum-energy configuration of $ in the xy plane is identical to the configuration of a 
point defect in the corresponding 0(2) vector model in d — 2. Thus, the structure factor is given by 

,L/2 ,L/2 

S(k) = $(k)-<f>{-k)=S {2 \k x ,k y ) / dze lk > z dz'e- lk ' z , (15) 

J-L/2 J-L/2 

where S^ 2 '{k x ,k y ) is the structure factor for the d — 2 point defect configuration. From Porod's law in d = 2 we 
have S^(k x ,ky) = /[V^ (k 2 + k 2 ) 2 }, where is a dimensionless constant, is the system area, and 9 is 
the angle between k and the orientation of the defect line. We now let L — > oo and calculate the structure factor per 
unit length of the defect, S scg (k) = linif^oo S(k)/L. By using the identity 
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L- 

we immediately obtain 



i f L/2 f L/2 
lim L 1 / dz exp(mz) / dz' exp(— iuz') = 27r<5(u) , (16) 

x J-L/2 J-L/2 



2 ttA^ 1 

V k 4 sir 



^c g (k) = T , - - , S(kcos9). (17) 



Eq. (17) can be used as the basis for the evaluation of the structure factor of an arbitrary vortex-loop configuration, 
provided that the local radii of curvature of the defect line are large compared to fc _1 at all points along the loop (i.e. 
the curvature affects the correlations of the order parameter only on length scales exceeding fc~ x ). We then reach the 
general result for the structure factor of a vortex loop in three spatial dimensions 

„ „ , 2nA^ 1 / , <5(fccos6»(s)) , x 

where s denotes the loop arc-length. For example, for a circular loop with radius R that has the normal to the loop 
plane oriented at an angle £ relative to k, we use the parameterization 9 = arccos[sin(£) sin(s/i?)] with s € [0, 2irR] 
and obtain 

„ „ . AnRA^ l 

V k | sin C| 

The fc -5 dependence in Eq. (|l^) agrees with the general Porod-law form, Eq. (0), with the exponent x = 2d — D = 5 
corresponding to the spatial dimensionality d — 3 and defect dimensionality D = 1. The result Eq. (|l^) is valid in 
the range R^ 1 <C k <C where £ is the core size of the disclination. Note that <S c irc reaches its minimal value, 
AttRA^ /Vk 5 , when k lies in the plane of the defect loop (i.e., £ = 7r/2), and diverges as the angle £ approaches or 
7r (i.e. k perpendicular to the loop plane). This reflects the fact that the length of the loop segments at which the 
angle 9 between k and the local vortex orientation is sufficiently close to tt/2 for the segment contribution Eq. jl7|) 



to be nonzero increases with decreasing £ (for ( = 0, we have 9 — it/ 2 all along the loop). Equation (19) ceases to be 
valid in the limit ( — > 0; for finite R, the regime of validity is restricted to angles |C| > n/(kR) |30| . 

It is worth stressing that there is no need for any large-distance cut-off in the loop calculations, as none of the 
integrals leading to Eq. ( |l9| ) diverges at small k. The factor of system volume V in Eqs. ( [T7jjT9| ) arises purely due to 
the chosen normalization of the structure factor. The presence for each loop segment of an antipodal loop segment at 
separation 2R is reflected in the structure factor result solely through the fact that the validity of Eq. ( JlQ| ) is restricted 
to k values larger than i2 _1 . 

Now consider the situation, typical during the phase-ordering process following a quench, in which the system 
contains many vortex loops of random orientations. In this situation, the properties of the system are globally isotropic 
and the structure factor S(k) depends only on the magnitude k of the scattering wave-vector. We may obtain S(k) 
by averaging the result in Eq. (|17| ) over all orientations of the defect segment: Sisotr(fc) = 1/2 Jq d9 sin 9 S seg (k, cos 9). 



By using the result (17) we obtain the structure factor per unit defect length in a globally isotropic system: 



7rA( 2 ) 1 

S isotI (k) = ^ r¥ . (20) 

So far, all arguments in this section have depended only on the scattering geometry (i.e. the shape and orientation of 
the defect line and the orientation of the wave- vector) , and were valid independent of the form of the order parameter. 
To obtain specific results for case of the 0(2) vector model in d = 3, it suffices to substitute the result A^ = Air 2 for 
the Porod amplitude of the point defect in the 0{2) vector model in d — 2 [Eq. (|lj) for d = N = 2] into the general 
expressions (|l7-20). For example, using A^ = Air 2 in Eq. (pFJ), we obtain the Porod amplitude in a phase-ordering 



0(2) vector system in d = 3 as A^ — Att 3 . This last result is in agreement with the corresponding result obtained in 
Ref. " 



IV. UNIAXIAL NEMATIC SYSTEMS 



A. The nematic order parameter 

For the case of nematic liquid crystals the local order-parameter field Q a ^(r,t) is a symmetric traceless rank-2 
tensor with Cartesian indices a, /3 = 1, 2, 3 (see, e.g., Ref. Q|). In the common case of the uniaxial nematic, the order 
parameter has the form 
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3 1 

Qa/3 — -^S\{u a up - -6 a p) , (21) 

where the order-parameter magnitude Si determines the strength of orientational ordering of the ncmatogcn molecules, 
and the director u gives the local value of the preferred orientation of the molecules. In the more general biaxial nematic 
case, the order parameter may be written as 

3 11 

Qap = -Sx(u a up - -8 af} ) + -S 2 (b a b f3 - v a v fj ) , (22) 

where ±u is the uniaxial director, ±b is the biaxial director, v = u x b, and the amplitudes S\ and S 2 determine, 
respectively, the strength of uniaxial and biaxial ordering. 

For the nematic, the real-space correlation function C(r,i) is defined as 

C(r) = — — f d 3 x Tr [Q(x) Q(x + r)] , (23) 

where Tr[A£?] = A a pBp a [cf. Eq. (j^)], and the normalization factor M ncm , which enforces C(r = 0) = 1, is given by 

d 3 xTr [Q(x) Q(x)] = ^(3S 2 + S 2 2 ) , (24) 

in which V is the volume of the system. For the nematic, Eq. (^) becomes 

S(k) = — L_ Tr[Q(k)Q(-k)]. (25) 

In this Paper, we are restricting our attention to unpolarized light scattering, in which case the structure factor ( pq ) 
is directly proportional to the measured scattered- light intensity [jy. 

When present, a given type of defect will adopt an equilibrium configuration of the nematic director that minimizes 
the nematic free energy. For future reference, we recall that the Frank free energy of a uniaxial nematic has (apart 
from surface terms) the general form [Q 

E = ^Jd 3 x [K n (V ■ u) 2 + K 22 (u • V x u) 2 + K 33 |u x (V x u)| 2 ] , (26) 

where the splay (Kn), twist (K22), and bend (-K33) elastic constants depend on the order parameter magnitude Si. 
Eq. (|2^) is often simplified by adopting the so-called one-constant approximation, Kn = K22 = = K, in which 
case the free energy (apart from surface terms) can be written as 

E = — J d 3 x(V lUj )(V lUj ) . (27) 



B. Hedgehog defects in uniaxial nematic systems 

A three-dimensional uniaxial nematic system admits topologically stable point defects (i.e., nematic hedgehogs) as 
well as line defects (i.e., nematic disclinations). (For a pedagogical discussion of the issue of topological stability, see, 
e.g., Ref. p]}].) The one-constant approximation to the Frank free energy, Eq. (pTj), admits (up to global rotations) 
two minimum-energy point-defect solutions with unit topological charge. In the radial hedgehog configuration, the 
director ±u points everywhere radially outwards from the center of the defect, and is given by u = r/r. The hyperbolic 
hedgehog configuration can be obtained from the radial hedgehog configuration by inverting one of the components 
of the director u, e.g. u = (—x/r,y/r,z/r). Note that the two configurations are topologically equivalent pi] ]. 

The calculation described below gives identical results for both the radial and the hyperbolic hedgehog; for simplicity, 
we shall work with the radial hedgehog configuration of the director, 

Q«/j(r) = |si (yy-^) • (28) 

The following calculation of S(k) closely mirrors the steps Eqs. (|)-((l|) in our O(N) vector model calculation. 
First, we observe that the Fourier transform of Q a p(r) is itself a symmetric, traceless, rank-2 tensor, and that due to 
symmetry it can depend only on the direction ±k. The general form of such a tensor is 



7 



3 / 1 
Q«/3(k) = -Si [k a kp - -5 a p k 2 ) g{k 2 



(29) 



where g(k 2 ) is a scalar function of the scalar k 2 . By contracting Eq. ( p9| ) with k a k@ and solving for g we can rewrite 
Eq. (ga) as 



k • r 



2 k 2 



„ik-r 



9Si 



,fc 2 



(2tt) j -<W 



<9A 2 



G(A 2 fc 



(30) 



In this expression the 5-function part corresponds to the forward-scattered beam, and will be dropped henceforth, as 
it does not influence the large-/c behavior. The function G is defined via 



G{X 2 k 2 ) 



drr e 



2 „iAkr 



(31) 



i.e., G is the three-dimensional Fourier transform of the potential r 2 , which can readily be shown to be given by 



C < XV) = % 



Evaluating the second derivative of G, and inserting it into Eq. (|3fj|), leads to 

Q a /3(k) 



9tt 2 Si / 1 2 

k a kp - -o^r k 



4fc 5 



->a/3 ' 



(32) 



(33) 



From Eq. ( |24| ) we see that the normalization factor M nem takes the value 3VS 2 /2. Combining Eqs. (p4|), ( |25| ) and 
(^) gives for the structure factor S(k) of a single nematic hedgehog defect in a system of volume V the result 



S(k) 



36tt 4 1 
V ¥ ' 



(34) 



Note tha t the Porod amplitude 367r 4 for the nematic hedgehog differs substantially from the Porod amplitude 127r 3 
(see Sec. Ill A) for the 0(3) monopole in d = 3. 



C. Disclination lines in uniaxial nematic systems 

We now calculate S(k) for the case of a disclination line defect in a uniaxial nematic system. In the topologically 
stable disclination configuration, the director rotates about the core of the disclination by ±180° (see, e.g., Ref. pl[). 
Disclinations with ±360° rotations are unstable towards the "escape in the third dimension" (3l) and do not have a 
singular core. We shall concentrate our attention on the topologically stable disclinations (appearing as thin threads 
in optical observations) as the ±360° disclinations (appearing as thick threads) do not give raise to a power-law 
contribution to S( k). 



In Section III B , we gave general geometric arguments that allowed us to express the structure factor of a linear 
defect in d = 3 in several specific configurations in terms of the Porod amplitude of the point defect in d = 2 
(obtained as the cross-section of the line defect in d = 3). Equations (p^-pc]) are valid regardless of the order parameter 
in question and may, thus, also be used for nematic systems. The remaining task therefore is to obtain an expression 
for A^> in the nematic case. 



As the director rotates about the core by 180°, it is clear that the approach of Sec. IV B , which is based on rotational 
symmetry, cannot be directly applied. It is readily seen, however, that the calculation of the Porod-law amplitude 
for the 180° nematic defect can be mapped on to the corresponding calculation for the 360° defect in the 0(2) vector 
model in d = 2. For the 180° uniaxial nematic defect, in the one-constant Frank free energy approximation, the 
order-parameter configuration is given (up to a global rotation) by Eq. ( plj ) with u(x) = ( cos |0(x), sin ^0(x), 0) , 
where <^(x) is the polar angle of the radius- vector x in the plane perpendicular to the disclination line. The correlation 
function Eq. (|23|) then becomes 
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C[H (r) = J d 2 x cos 2 Q^(x) - ^(x + r)) + c a (35a) 

'• iS,/2h ^ilcos^x) -0(x + r)) + c 2 , (35b) 

where c\ and C2 are numerical constants. On the other hand, for the 360° defect in the 0(2) model in d = 2, the 
order-parameter configuration is given by <fr(x) = ( cos </>(x), sin </>(x)) , and the corresponding correlation function, 
Eq- is given by 

Cill (r) = ~jTj~ / d 2 x cos[0(x) - 0(x + r)] . (36) 

After Fourier-transforming Eqs. ( |35b) ) and ([56j) and omitting the 5-function term arising from the constant ki in 
Eq. (^), we obtain that the structure factors of the two defects are related simply by 

5l80 o (^) = ^77 ' "^360° (^) ' (37) 

^"nem 

Using the result S^ QO (k) = Att 2 ^ 4 /V^ [i.e., Eq. (Q for d = 2, with the area of the system denoted by U (2) ), and 

the normalization factors M ( 2 ) — an d Mncm — V^3Sf/2, we finally obtain the structure factor s[^ 0o (k) for 
the 180° point defect in a two-dimensional nematic: 

By using this result (i.e., = 3vr 2 ) in Eq. ©, we obtain the structure factor of a nematic disclination loop of 
radius R 

„ 127T 3 i? 1 

gcircfc = 39 

K fc° | sm(| 

Likewise, by using Eq. ( pp[ ) we obtain the orientation-averaged contribution to the structure factor from a segment of 
length L of a nematic disclination 

^) = |3- 3 ^- (40) 

It follows that the structure factor for a globally isotropic system containing disclination lines is given by S(k) = 
/9 37r 3 fc~ 5 , with the defect density p having the meaning of the total disclination length per unit volume of the system. 



D. Ring defects in uniaxial nematic systems 



Topologically stable disclinations lines in a uniaxial nematic cannot terminate in the bulk — they must form closed 
loops, or else extend to the system boundary. Depending on the details of the director configuration, the disclination 
loop can carry any (integral) monopole charge [ p2[ . In this section, we discuss typical examples of disclination loops 
with zero and nonzero monopole charges and their contributions to the structure factor. 

Consider first the case of a "twist disclination" loop, which carries zero monopole charge. The cross-section of 
the director configuration around a twist disclination is shown in Fig. |l|a. The axis of rotation of the director is 
perpendicular to the direction of the disclination line. The director configuration obtained by forming a closed loop 
of a disclination that locally has the "twist" structure of Fig. |l|a is illustrated in Fig. [j]b. This director configuration 
is homogeneous at large distances from the loop, and consequently does not result in a power-law contribution to the 
structure factor for k -€i For length scales smaller than R (i.e., for k » R^ 1 ), the twist disclination loop is 

characterized by the structure factor given in Eq. (39). 

Next, consider the "wedge disclination" loop. In this configuration, the director rotates about an axis parallel to 
the disclination lin e (see Fig. |2|a). The resulting director configuration (see Fig. ||b) has the structure of a radial 
hedgehog (see Sec. IV B| ) outside of the disclination loop. This results in a power-law contribution to the structure 
factor characterizing the nematic hedgehog in d = 3, Eq. (|34|), valid for k < R~ l , where R is the ring radius. On the 
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(a) (b) 

FIG. 1. Order parameter configuration around a twist-type ring defect in a uniaxial nematic system, (a) Cross-section of 
the director configuration around a disclination segment running perpendicular to the page with the core passing through the 
point "x" . The "nail head" marks the end of the director that lies in front of the page. The director rotates by 180° around the 
core of the disclination; the rotation axis is vertical, i.e., perpendicular to the orientation of the disclination segment, (b) The 
director configuration in the plane of the twist disclination loop. The director outside the disclination loop is rotated by 90° 
with respect to the director inside the loop; the rotation axis is everywhere perpendicular to the loop plane. The cross-section 
shown in (a) corresponds to the configuration around the disclination segment marked by "x" in (b). 




(a) (b) 

FIG. 2. Order-parameter configuration around a wedge-type ring defect in a uniaxial nematic system, (a) Cross-section of 
the director configuration around a disclination segment running perpendicular to the page with the core passing through the 
point "x" . The director rotates by 180° around the core of the disclination, with the rotation axis parallel to the orientation 
of the disclination segment, (b) The director configuration of the wedge disclination loop. The director is perpendicular to 
the loop plane everywhere inside the loop. Outside of the loop (in any plane), the director adopts the radially symmetrical 
configuration of a nematic hedgehog. The cross-section shown in (a) corresponds to the configuration around the disclination 
segment marked by "x" in (b). 
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other hand, for k 3> R l , we have the usual disclination-loop structure factor, Eq. (|39|). To summarize, the structure 
factor of a wedge disclination loop is given by 



S(k) 



36^ 1 

V W 
12tt 3 R 1 
V k 5 | sin C| ' 



k < R- 1 ; 

r- 1 <fc<r 1 



(41) 



Therefore the large-fc region of the structure factor measured in a system containing a wedge disclination loop of 
radius size R is expected to exhibit a crossover between two power-law regimes with exponents 6 and 5 (Fig. |3|). The 
crossover point, as defined in Fig. |^ , is predicted to occur at the value of k c = 3tt(/R. Note that measuring the 
location k c of the crossover allows one to obtain information about the ring-defect size R from un-normalized data 
for S(k); this should be contrasted to the use of the Porod law based on Eq. (40), where the total string length is 
extracted from the (rarely experimentally available) normalized structure factor. 
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k ( (o.m ) 



FIG. 3. Structure factor 5*(k) of the ring defect configuration in Fig. |^. The scattering wave-vector k is assumed to be 
oriented at an angle C 7^ with the respect to the ring axis. The two Porod regimes, arising for length scales larger or smaller 
than the ring radius R, are shown as the dot-dashed line [S(k) = 36tt 4 / (V k 6 ')] and the dotted line \S(k) = 12ir s R/(Vk 5 \ sinCD], 
respectively. The solid curve schematically shows the crossover between the two power-law regimes. The crossover wave-vector 
k c , as defined in the picture, is given by k c = 37ri? _1 | sin£|. Thermal fluctuations can further modify the form of S(k) sh own i n 
the figure, if the distance to the nearest other defect or the system boundary is large; these effects are discussed in Sec. [IV 3 . 

Even though the wedge disclination loop is a topologically admissible configuration, the (energetics-dependent) 
question of its occurrence in real nematic systems has recently attracted some interest. The energetic stability of the 
wedge loop configuration was investigated theoretically in Refs. and 34 1. There, it was found that provided that 
certain restrictions on the elastic constants in the nematic free energy are satisfied, there does indeed exist a non-zero 
equilibrium radius of the loop: a ring of a larger (smaller) radius will tend to shrink (expand). Depending on the values 
of the elastic constants, Ref. j34|] predicts i? cq in the range 10£ — 10 4 £, where £ is the coherence length of the nematic 
order parameter (i.e., £ is of the order of 10 — 100 A). Note that in a nematic with a positive dielectric or diamagnctic 
susceptibility anisotropy, the loop can be made unstable towards expansion to a larger radius by applying an electric 



or magnetic field along the axis of the defect ring 
were recently observed in numerical simulations 



51. Wedge disclination loops with a non-zero equilibrium radius 
Experimental observations of such a configuration, however, are 
currently not available. The form of the structure factor predicted in Eq. ([ll]), with its characteristic crossover, may 
be used as an experimental signature for the observation of the wedge ring configuration |52] . 
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E. Thermal fluctuations of director orientation 



Up to this point, we have ignored thermal fluctuations — we have effectively assumed that the structure factor 
was measured at zero temperature. Nematic liquid crystalline phases, however, typically occur and are studied in the 
room-temperature range, and thermal effects can play an important role the scattering of light. Due to the presence of 
low-energy fluctuations in the orientation of the nematic director, the system exhibits enhanced turbidity throughout 
the whole range of temperatures where the nematic phase is stable. These transverse fluctuations of the director 
result in a power-law contribution to the structure factor S(k), even when no defects are present in the system. 

We now very briefly review the theory of scattering by thermal fluctuations in bulk uniaxial nematic liquid crystals 
|fl, 22|. Consider a well-aligned region of the sample with the average director Uo parallel to the z-axis. The local 
va ue of the director can be written (to first order in uj_) as u(r) = Uo + uj_(r), where the transverse fluctuation 
uj_(r), by definition, satisfies Uj_(r) _L Uo and we assume that |uj_(r)| -C 1. The nematic structure factor, Eq. (pq), 
then becomes [upon dropping terms proportional to 5(h)] 

S(k) = A u±(k ). ux (_ k ). (42) 

Note that the quantity uj_(k) • u^(— k) must now be regarded as a thermal average. It's value is readily obtained 
from the equipartition theorem. In the "one-constant" approximation, using u(r) = Uo +u^(r) in Eq. ( p7| ) results in 
the free energy 

1 f d 3 k , 
E=-Kj—k^)\\ (43) 

Furthermore, it is possible to decompose uj_(k) into two orthogonal modes that decouple Q); by equipartition, we 
then have |u^(k)| 2 = k^TV. Equation (|4|) thus becomes 

*»>-il&- < 44 » 

Taking a typical value K = 10~ 6 dynes for the elastic constant and k-gT ~ 0.5 • 10 -13 erg for the thermal energy at 
room temperatures, one obtains the estimate 

S T (k)^^. (45) 

Note that in the light-scattering literature, the structure factor is often defined without the normalization factor Af nem 
in Eq. (^5|), and then has the dimension of volume squared, rather than of volume as in our case. 

It should be stressed that Eq. ([S|) was derived from a free energy appropriate for small fluctuations in a well-aligned 
part of the nematic sample. Equation (Q) can therefore be used to describe scattering from a system with a spatially 
varying director only provided that there are no appreciable variations in the static director orientation on the scale 
offc- 1 . 

Even in strongly inhomogeneous samples containing topological defects, thermal fluctuations will dominate the 
scattering intensity provided either that the defect density is sufficiently low or that k is sufficiently high. To illustrate 
this point, we now consider light scattering from a nematic droplet of radius R with normal (homeotropic) boundary 
conditions on the droplet surf ace. I n the minimum-energy configuration, the droplet contains a radial hedgehog in the 



center. (As discussed in Sec. IV D, the radial hedgehog is, in fact, expected to take the form of a wedge disclination 
loop of small radius r; in the present paragraph, we restrict our attention to length scales exceeding r and do not 
consider the detailed structure of the hedgehog.) Consequently, for a given scattering wave- vector k 3> 2ir/R, the 



central region of the droplet will contribute to the structure factor according to Eq. (34), Sa c f(k) = V~ lr S6ii A k~ & , 
where V — AttR 3 /3. In the region 2ir/k < r < R, the static order-parameter configuration does not contain variations 
that contribute substantially to S(k); however, thermal fluctuations on the scale 2ir/k still exist. The contribution 
to S(k) from these is given by Eq. ((45|). The resulting ratio j(k) of thermal and static scattering intensities, i.e. 
7(fe) = ST{k)/ Sdef(k) = 0.04(i?fc) 3 fc, where k is given in inverse A, can very in a wide range. Clearly the static 
contribution Sdci(k) dominates when k is sufficiently small (i.e. near the forward-scattered beam). Consider, on the 
other hand, scattering at the experimentally accessible value of k = 10 _5 A _1 . For a droplet of radius R = 1 cm, we 
obtain 7 ~ 400; for a droplet of smaller radius R — 100 fim, we obtain 7 ~ 0.0004. We see that either the thermal 
or the static contribution to the scattering intensity can completely dominate for R and k in the experimentally 
reasonable range. (For the sake of completeness we note that we have not taken into account scattering from the 
droplet surface in the arguments given above.) 
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F. Porod tails in the phase-ordering kinetics of uniaxial nematics 



We are now ready to discuss in detail the short-distance form of the structure factor arising during the phase- 
ordering process in uniaxial nematic liquid crystals, and to relate our conclusions to the experimental results reported 
in Refs. §,§. 

We first briefly recall some general features of phase-ordering kinetics following a quench from a disordered to an 
ordered phase [ p7y3q| . Phase-ordering kinetics initially attracted theoretical attention due to the property of dynamical 
scaling of the order-parameter correlations at late times after the quench. In its simplest form, the corresponding 
dynamical scaling hypothesis states that all time-dependent length scales in the system have the same asymptotic 
time dependence (and, consequently, one can define a single "characteristic" length scale). In most systems with 
non-conserved order-parameter dynamics, this time-dependence is given by the power law L(t) oc t 1 / 2 . In systems 
containing topological defects, the characteristic length scale L(t) can usually be extracted as the typical defect-defect 
separation at time t. 

The process of phase ordering has been successfully studied experimentally in systems described by a scalar order 
parameter (e.g., in binary alloys — see Ref. [(39[ ). However, analogous experiments were found to be difficult to 
perform in most systems with continuous symmetries (such as ferromagnets or liquid He 4 ). In the recent years, a 
series of experiments f|,|9|,[5],[| have demonstrated that nematic liquid crystals provide a system in which phase ordering 
(following a quench from the isotropic to the nematic phase) is readily accessible to experimental investigation. In 
Refs. the phase-ordering process was studied by measuring the time-dependent nematic structure factor S(k, t). 
The present section includes a discussion of the large-fc region of the structure factor measured in these experiments; 
however, we also discu ss effe ct s tha t may be observed under different experimental conditions. 

As we saw in Sees. ( IVB )-( [VD|) , a uniaxial nematic liquid crystalline system can contain, simultaneously, point 



defects (hedgehogs) and line defects (disclinations) |4(J . If the average separation between hedgehogs is given by Lhedg 
and that between disclinations by Ldisci, a unit volume of the system contains phedg = ^hedg hedgehogs and a length 
Pdisci = disci °f disclinations. Now consider scattering at a wave-vector k such that fc _1 <C Ldisci and fc _1 <C Lhedg- 
The total structure factor per unit volume due to hedgehogs is then given by phedg 367r 4 /fc 6 [see Eq. [m| , and that due 
to disclinations by /?disci 37r 3 /k 5 [seeEq. ©]. (Here we have assumed that the disclination lines are predominantly of 
the twist type, as is expected from energetic considerations.) In the parts of the system that are not within distance 
k^ 1 from the nearest disclination or hedgehog (which, by assumption, is most of the volume of the system), the static 
director configuration is effectively homogeneous on scales of order of ft -1 . Consequently, Eq. ( J44| ) for the strength 
of scattering due to thermal fluctuations applies, and per unit volume, we have St(&) ~ 8A/k. (Notice that the 
length scale occurring in the expression for Sr{k) is microscopic, as opposed to the defect separations occurring in 
the hedgehog and disclination contributions.) The total structure factor is therefore given by 

3671-4 37r3 8 ^ 

S(k) = p hodg — + Pdisci -p- + -p- • (46) 

First consider a system in which disclinations dominate over hedgehogs. We now estimate, for a given density of 
disclinations pdisci (and phodg = 0), the range of scattering wave- vectors k in which the the generalized Porod law 
behavior, S{k) oc A: -5 , can be observed. The lower limit on the range of k is given by the inverse average separation 
of defects: k\ = L~^ scl — p\ iscV The upper limit fc u is given by the k value for which ST(k) = 8Ak~ 2 exceeds the 

Porod law contribution, Sdisci — Pdisci37r 3 fc~ 5 . This gives k u — (Sir 3 /8A) 1 / 3 L^ 2 ^ 3 . We see that although both k\ and 
k u decrease when the separation of defects Ldisci increases, the ratio k u /k\ increases as (i disc i/0.09l) 1/3 . Analo gous 
estimates for observing the S(k) oc fc~ 6 behavior in a system dominated by hedgehogs yield fci = i^edg = /°hcd g > 

fc u = (367T 4 /8A) 1 / 4 L^ c 3 d / 'g , and k u /k\ — (Li lc d g /0.002A) 1 / 4 . Consequently, the Porod form of the structure factor is 
predicted to extend over one to three decades in k for the experimentally accessible values of defect separation of the 
order of 0.1 /im to 1cm, for either a disclination-dominated or a hedgehog-dominated system jyj. Crossover effects 
between fc~ 6 and fc~ 5 scaling in situations when both hedgehogs and disclinations are present may further diminish 
the range in which these power laws can be clearly observed. (A specific example of such behavior will be seen below.) 

It should be noted that measuring the location of the crossover wave-vector fc u permits the efficient estimation of 
the density of defects in the system. Specifically, if one measures the location of the crossover k u in a disclination- 
dominated system, the total disclination length per unit volume (expressed in units of A ) is given by (2.3/fc u ) 3 / 2 . 
In a hedgehog-dominated system, the number of hedgehogs per unit volume (expressed in units of A ) is given by 
(4.6/fc u ) 4 / 3 . Although it is also possible to estimate the defect density from the position fci of the crossover at low 
wave- vectors, it should be noted that this crossover is less steep than the crossover at fc u , and its precise character 
depends on the nature of the correlations amongst the defects. 
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We now relate our results to the measurements of S(k) in a uniaxial nematic undergoing phase ordering in the 
experiments of Wong et al. ||||. The range of scattering wave-vectors investigated in Refs. was fc = 1000cm -1 - 
10000cm -1 , and the typical defect separation in the system at intermediate and late times after the quench was of 
the order of 10/itm (see Fig. 1 in Ref. Direct optical observations in this and related |||| phase-ordering systems 
indicate that disclinations dominate over hedgehogs; the exact relative proportion of their densities is, however, 
difficult to measure. 

10" 5 | ■ ■ ■ ^-| 




FIG. 4. Theoretically predicted short-distance structure factor S(k) for a 3-dimensional uniaxial nematic system undergoing 
phase ordering. The lower solid curve shows S(k) for disclination density pdisci = W 6 cm~ 2 and hedgehog density zero. The 
upper solid curve corresponds to disclination density pdisci = 10 6 cm -2 and hedgehog density phedg = 10 s cm -3 . The long-dashed, 
dashed, and dot-dashed lines show the contributions due to disclinations, hedgehogs, and thermal fluctuations, respectively. 



In Fig. [|, we plot the structure factor S(k) predicted by Eq. ( fig ) for two situations in which the average separation 
between defects is of the order of 10 /im. The lower solid curve shows S(k) for pdisci = 10 6 cm -2 and phedg = (i.e., a 
configuration dominated by disclinations with Ldisci = 10 /im). The upper solid curve corresponds to pdisci = 10 6 cm -2 
and Phedg = 10 8 cm -3 (i.e., the system contains, in addition, hedgehogs with average spacing Lhodg = 22 /im). We 
restrict the plot to the region k > (10 /im) -1 = 10 3 cm -1 , where Eq. ( [iq ) is expected to be applicable. The broken 
lines in Fig. || show the three individual contributions to S(k) arising from hedgehogs, disclinations, and thermal 
fluctuations. 

It is seen that for both phedg = and phedg = 10 8 cm -3 , the crossover to the thermal-fluctuation-dominated regime 
[S(k] oc fc -2 ] occurs at approximately fc = 50000 cm -1 , i.e., well beyond the range of k measured in Refs. |||6| and at 
the limit of the range of k accessible, in principle, with visible light. Indeed, the structure-factor data in Refs. ^|J|] 
do not show any sign of such a crossover. It should be noted, however, that if the scattering experiments would have 
been performed at a (still later) time when the average defect separation would have reached 300 /im, the crossover 
value would have decreased to k ~ 5000 cm -1 and the sharp crossover to the thermal regime would have been clearly 
visible with the experimental setup of Refs. Hp] provided that the scattering probe were sufficiently sensitive. 

In the case phedg = 0, the structure factor in Fig. || exhibits the Porod behavior S(k) oc fc -5 , characteristic of 
disclinations, over a range of 1.5 decades (fc = 1000cm -1 to fc = 50000cm -1 ). In the case phedg = 10 8 cm -3 , a 
crossover between S(k) oc fc -6 and S(k) oc fc -5 is spread between fc = 1000 cm -1 and 10000 cm -1 , and leaves only a 
very narrow range of fc for which the S(k) oc fc -5 scaling behavior may be observed. The fact that hedgehogs result 
in a substantial modification of the structure factor even when iphedg/pdisci = 0.1 is rather surprising, and is due to 
the large ratio, Ahodg/Aiisci = 127r, of the Porod amplitudes for hedgehogs and disclinations. 

The structure-factor data reported in Refs. J|,U (for three-dimensional samples) were fit in these references to 
a power-law with exponent 6 ± 0.3. Some of these data (specifically Fig. 7 in Ref. were later re-analyzed in 
Ref. 0, with the conclusion that the asymptotic slope was closer to 5, but was approached from above , i.e., through 
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effective exponents lying between 5 and 6. Our results show (see the upper solid line in Fig. ||) that such behavior 
is expected to arise if a sufficient number of hedgehogs is present in the system. The relative density of disclinations 
and hedgehogs was not examined in the experiments of Refs. Other experiments investigating phase ordering in 

nematics, however, indicate that hedgehogs are rare compared to disclinations. (For the sake of completeness we note 
that it is possible to prepare nematic systems dominated by hedgehog defects: thus, in Ref. [Q, a liquid crystalline 
system was quenched from the isotropic to the isotropic-nematic biphasic region, and a large number of hedgehogs 
formed upon coalescence of the nematic droplets.) Specifically, it was found in Ref. |)| that hedgehogs occurred in 
significant numbers only within a specific intermediate range of times after the quench, and that the ratio Lphedg/pdisci 
within this time range was of the order of 10-100. The case pdisci = 10 6 cm -2 , phedg = 10 8 cm~ 3 , chosen above to 
obtain the upper solid curve in Fig. |J, falls within this range. Consequently, it cannot be ruled out that the observed 
approach to S(k) cx fc~ 5 from above in Fig. 7 of Ref. || is due to the presence of hedgehogs, even though if the 
relative proportion of populations of hedgehogs and disclinations in the experiments of Refs. |5|(| is similar to that 
found in the experiments of Refs. || (which was performed on a different nematic material), such a possibility should 
be considered unlikely. A more l ikely possibility is that a portion of the length of the disclinations in the system is 



of the wedge type (see Sec. |VD|). Any curved configuration of a wedge disclination will lead to an order-parameter 



configuration in the vicinity of the disclination that resembles a part of a (radial or hyperbolic) he dgeho g 

m- (The 



special case of a disclination loop that is wedge-like everywhere along its length was discussed in Sec. IV D ; in this case, 
a full hedgehog configuration is obtained outside of the loop.) Although the twist configuration of a disclination line 
is energetically preferable to the wedge configuration for usual values of the nematic elastic constants |Q , wedge- type 
segments of disclination loops in a system undergoing phase ordering may be generated dynamically, and lead to a 
k~ & contribution to the nematic structure factor. 

We complete this Section by a brief discussion of the dynamical scaling of the structure factor. Three different 
length scales — the disclination separation I/disci, the hedgehog separation Lhcdg, and the microscopic length scale 
8 A associated with thermal fluctuations — appear in the large-fc form of the nematic structure factor, Eq. (pi|). Is 
this inconsistent with the property of dynamical scaling of the structure factor, according to which there exists a 
(time-dependent) length scale L(t) such that k 3 S(k) depends only on the dimensionless scaling variable y = kL(t)l 
Let us choose the disclination separation Ldisci as the characteristic length L(t). From Eq. (p6[), we then have 



fc 3 ^) = 3^J_ + ^[W!) 

ihedg(i) 



^disclitj 



Clearly, k 3 S(k) does not in general satisfy the scaling form /(fcLdisci(^))- At asymptotically late times, however, 
^disci 3> 8A and, consequently, the term in Eq. ([!?]) associated with thermal fluctuations becomes negligible. Ex- 
perimental evidence |9) indicates that at late times after the quench, the hedgehog separation Lhcd g (£) grows faster 
than the disclination separation Ldisci(i)- Consequently, the second term in Eq. ( |47j ) (associated with hedgehogs) also 
becomes negligible, asymptotically. Our analysis is therefore consistent with the dynamical scaling of the structure 
factor at asymptotically late times. As we saw earlier in this Section, however, the scattering contributions from 
thermal fluctuations and from nematic hedgehogs can lead to significant transient effects under typical experimental 
conditions. 



V. BIAXIAL NEMATIC SYSTEMS 



Thus far, our discussion of defects in nematics has concentrated on the case of uniaxial systems. It is straightforward 
to generalize this discussion to the case of biaxial systems. As there are no topologically stable point defects in a 
three-dimensional biaxial nematic p5| , we need only consider line defects. Biaxial nematics admits four topologically 
distinct classes of line defects ji5| , which are disclination lines distinguished by the angles of rotation of the uniaxial 
director u and the biaxial director b [defined in Eq. (p2|)j around the defect core. In the C x class of disclinations 
u rotates by ±180° and b does not rotate; in C y disclinations u does not rotate and b rotates by ±180°; in C z 
disclinations both u and b rotate by ±180°; finally, in Co disclinations either u or b (or both) rotate by 360°. 
These four distinct disclination types were observed experimentally in a thermotropic nematic polymer |jl9[] , and their 
properties where found to be in agreement with the predictions of the topological classification scheme. 

In the minimum-energy configuration [with free energy given by Eq. (^?J for the C x , C y , and C z defects, the 180° 
rotations of the u and b director are uniform f46f| . The configuration of a C x disclination is correspondingly given 
(up to a global rotation) by Eq. ( p2| ) with 

u(x) = ( cos ~<£(x), sin ^(x), 0) , (48a) 
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b(x) = (0,0,1), 

1 7T 1 

v(x) = ( cos (-</>(x) + -) , sin (-< 



4 y 



(48b) 
(48c) 



where </>(x) is the polar angle in the plane perpendicular to the disclination line. The cor rela tion function (23) for this 
configuration can be expressed (up to an additive constant) as the correlation function (35t) for the uniaxial nematic 
disclination, multiplied by the biaxiality-strength-dependent factor i?(35i +S , 2) 2 /4, where R = 35 2 /(35 2 + 5|) is the 
ratio of the uniaxial (52 = 0) and biaxial (S2 ^ 0) normalization factors M neia . The corresponding factors for the case 
of the C y and C z disclinations are likewise readily evaluated. By using the result (40) for the uniaxial disclination, 
we obtain the structure factors of the biaxial disclinations C x 



C y and C z 



S(fc)c x 
S(k) Cy 



9tt 2 L x 5?(35i + S 2 ) 2 1 



8 V 
9tt 2 L v 



3Sf-\ 

o2 c2 
°1 °2 



si 



k 5 



2 V 35 2 + 
9tt 2 L z 5 2 (35 



Stk* 



S2) 2 1 



V 35 2 + si 



k 5 



(49a) 
(49b) 
(49c) 



It should be noted that in the case 35i = 52, the amplitude in the Porod law for the C z d ef ect is zero, whereas the 
amplitudes for the C x and C y defects are non-zero and equal to each other [see Eqs. (49a)-(49c)]. This reflects the 
fact that for 35i = 52, the order parameter given by Eq. (p2) describes a uniaxial discotic phase |47[ with uniaxial 
axis v = u x b. Correspondingly, the C z configuration in this case does not represent a defect (as v does not rotate), 
whereas the C x and C v configurations are equivalent. 



It remains to consider the Co disclination. For the case 5i > 52 (i.e., needle- like ordering 47 1), the minimum-energy 
configuration of type Co is given by Eq. ( p2[ ) with 

u(x) = (0,0, 1), (50a) 
b(x) = (cos0(x),sin0(x),O), (50b) 
v(x) = ( cos (0(x) + 7r/4) , sin (<j>(x) + tt/4) , 0) , (50c) 

where <^(x) is the polar angle in the plane perpendicular to the disclination line. The correlation function ( |23[ ) for 
the point defect in the planar cross-section through this disclination can then be expressed as i?(252/35i) 2 , By using 
Eq. ( pfl| ) one obtains the structure factor per unit length of the Cq disclination: 



12vr 



2 L 



Si 



S( * k)c ° ~ ""' V 35 2 + S 2 2 k 5 



(51) 



For the case Si < S2 (i.e., discotic ordering p7[]) the v director (and not the u director, in contrast with the needle-like 
case) corresponds to the eigenvalue of the order-parameter tensor Q a p largest in absolute value, and the lowest energy 
configuration that can be taken by the Co disclination has directors given by 



u(x) 



cos </>(x) , sin 0(x) , 0) , 

7T\ 



b(x) = ( cos (<Xx) + -) , sin (0(x) + -) , 0) , 
v(x) = (0,0,1). 



The corresponding structure factor is 



8L 2 (35 2 - 5 2 ) 2 1 



35 2 



Si 



fc 5 



(52a) 
(52b) 
(52c) 

(52d) 



In general, t he la rge- k structure factor of a biaxial nematic system is given by the sum of the expressions (49a-49t) 
and (H) or ( fgg ). 

Finally, we note that in a biaxial nematic film (i.e., a system having two spatial dimensions, but the full 3x3 
tensorial nematic parameter), C x , C y , C z and Co are point-like defects (the dynamics of which was studied in po|). 
The Porod-law amplitudes for these defects are obtained from the amplitudes in the results (49a)-(49c), ( |5l| ) and 
(^32d|) by dividing by w [see Eq. (p0|)], and the corresponding Porod-law exponents take the value 4 instead of 5. 
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VI. CORRECTIONS TO POROD'S LAW 



In this section we briefly discuss several effects — inter-defect interactions, inequality among the three elastic con- 
stants in the nematic free energy, curvature of linear defects, and the presence of the defect core — that were not fully 
taken into account during our calculations in the previous sections and may, under certain circumstances, lead to 
modifications of our results. 

An interesting question, already alluded to in Sec. ||, is whether the order-parameter configuration in the vicinity 
of the defect core is or is not affected by the interactions with the other defects present in the system. [In order to 
avoid confusion with the core region, where the order parameter magnitude is reduced, we shall refer to the region 
close to, but outside, the core, as the "central" region of the defect.] It is important to distinguish between static 
and dynamic effects. In the case of the 0{2) vector model with the standard gradient free energy, given by Eq. (|^), 
the minimum-energy configuration of a collection of defects with fixed locations is obtained simply by taking the 
superposition of the angles characterizing the order parameter m around isolated defects; their central regions are 
therefore unaffected. A deformation does occur, however, once the defects are allowed to move flq| ; for sufficiently 
slow defect velocities, this effect may be neglected. 

The situation is different in the 0(3) vector-model case. Here, Ostlund |^6| showed that in the minimum-energy 
configuration of a monopole-antimonopole pair separated by a fixed distance, the gradient free energy is concentrated 
along a string connecting the two point defects. We are not aware of any investigation of how this result might be 
modified in a dynamical situation. Ostlund's result does suggest, however, that the order-parameter configuration 
of monopoles and antimonopoles in a phase-ordering 0(3) vector system in d — 3 may substantially differ from the 



radially symmetric configuration assumed during the calculation of the corresponding Porod law in Sec. [II A of the 
present Paper. 

Consider, then, an arbitrary point defect in the O(N) vector model in d = N dimensions. We denote the corre- 
sponding order-parameter configuration by 3?(r), where the radius- vector r originates at the core of the defect. As 
already indicated in Sec. [n|, the value £ = 2d of the Porod exponent for point defects is suggested by simple dimen- 
sional analysis, independently of the exact form of 3>(r). This argument can be made more precise if we consider 
configurations in which 3>(r) depends only on the direction f, and not the magnitude r, of the radius- vector r. We 

may write <fr(k) = J d d r e ik ' r <&(r) = J d d r e lfck ' r $(r); by assumption, $(fc _1 r) = «J?(r), which, combined with the 
substitution r = fc _1 y, yields «fr(k) = k~ d J d d y e 4k ' y 3>(y) = k~ d <fr(k). Consequently, the structure factor Eq. (^) is 
of the form S{k) = A{k)k- 2d , where the Porod amplitude A(k) = M"^*^) ■ *(-k) now depends on the orienta- 
tion of k. The Porod exponent £ = 2d, expressing how S(k) depends on the magnitude of k, is, however, unchanged 
compared to the radially symmetric case. As for the value of ^4(k), we may make the following argument using the 
identity relating the standard form of the gradient free energy E to the structure factor S(k): 

E=^J dtxV&iWV&fc) = lj d d kk>S(k) = |(A(k)) ang | c i dkk^ d . (53) 

Here (^4(k)) ang denotes the angular average of A(k) over all orientations k, k is a positive elasticity constant, £ is the 
core size of the defect and D is the large-distance cutoff necessary in the calculation of the free energy. Equation ( p33] ) 
implies that (A(k)) ang reaches its minimum for the order-parameter configuration $(f ) that minimizes the free energy 
E |ft9j . In the case of the O(N ) mon opole with charge +1, the configuration minimizing E is precisely the radial 



configuration used by us in Sec. [II A, and we conclude that our result Eq. ( |14| ) provides a lower bound for the value 
of the Porod amplitude in a macroscopically isotropic system f30|| . 

Note, however, that the property of (^4(k)) ang formulated above applies equally well for defects that do not possess 
radial symmetry in the minimum-energy configuration, such as monopoles with higher-integer charges in the O(N) 
vector model, or strength- 1/2 disclinations in nematic liquid crystals [in the latter case, E in Eq. ( p3|) is replaced by 
the Frank free energy in the one-constant appro ximation, Eq. (E7I)] . We thus reach the conclusion that in nematics for 



which the three elastic constants (see Sec. [V A ) are not equal to each other and consequently, the configuration of the 



nematic director around a disclination differs from the configuration minimizing Eq. (ETU, the angle-av eraged stru cture 



factor in a phase-ordering system will have Porod amplitudes larger than those derived by us in Sees. IVB-IVF, even 
when the central regions of the defects are not deformed by inter-defect inter actions. 



While deriving the results for linear defects [0(3) vortex lines in Sec. [II B and nematic disclinations in Sees. IV C 



and IVD], we assumed that the curvature 1/R of the defect lines was negligible (i.e. Rk 3> 1). In general, we can 
expect corrections to the basic Porod law k~ 5 that are suppressed by factors of Rk, thus resulting in contributions to 
S(k) decaying as hig her power laws (fc -1 ' with r\ > 5). The case of curved wedge- type disclinations (discussed near 



the end of Sec. 1VC), for which r\ = 6, serves as an example. 
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Similarly, we can expect higher-power terms to arise from the presence of the core region of the defects. As the 
crossing of the boundary between the core region and the outside region of a defect is associated with a change in 
the magnitude of the order parameter, we may obtain an upper-bound estimate of the corresponding contribution 
Score (fc) to the structure factor by regarding the core boundary as a domain wall in a scalar system, leading (in a 
3-dimensional system) to S COIC (k) ~ afc~ 4 , where a is the domain wall area. For a monopole a ~ £ 2 , where £ is the 
core size, and the ratio of S COTe (k) to the standard monopole contribution k~ 6 is of the order of (£fc) -2 . For a string 
defect of length L, we have a ~ £L, and the ratio of 5 l C orc(fc) to Lk~ 5 is of the order of £fc. For typical core sizes £ 
of the order of 10A and for scattering wave-vector values k < 10~ i A accessible with visible light, we have £/c > 1000, 
and the contributions to S(k) from the defect core are expected to be negligible for both monopoles and strings. 



VII. CONCLUSIONS 



In this Paper we have presented a detailed discussion of the influence of topological defects on the large- wavevector 
behavior of the structure factor S(k) in nematic liquid crystals. The presence of topological defects leads to power- 
law contributions to S(k) of the Porod form pAk~^, where p is the number density of a given type of defect, A is a 
dimensionless amplitude, and £ is an integer-valued exponent. We have computed the values of the Porod exponents 
and amplitudes for the various types of topological defects present in uniaxial and biaxial nematics, and have discussed 
the competition between contributions to 5(k) due to defects and to thermal fluctuations. 

Our main results are summarized in Tables [| (for nematic films) and |fi] (for bulk nematic systems). To obtain the 
short-distance structure factor of a nematic system containing topological defects, the expressions given in Tables || 
or H should be multiplied by the number densities of the corresponding defects and then added. Here, the defect 
number density is defined as the number of point defects, resp. the total length of line defects, per unit area (in 
a two-dimensional system) or unit volume (in a three-dimensional system). In addition, the total structure factor 
contains a power-law contribution due to transverse thermal fluctuations of the nematic director, this contribution 
being given (in a bulk system) approximately by 8A/k 2 . The resulting form of the structure factor S(k) is valid for 
k ranging from the inverse typical separation of defects to the inverse defect-core size. 

To avoid confusion, we remind the reader that our results were calculated for a structure factor S(k) defined so 
that the real-space correlation function C(r) is normalized to unity at r = (and, consequently, our structure factor 
has the dimension of volume). Thus, our results differ from the un-normalized structure factor (with dimension of 
volume squared) by the normalization factor M n0 rm = (V/2)(SSf + Sf) (where V is the system volume or area, Si 
is the (uniaxial) order-parameter magnitude, and S2 is the strength of biaxial ordering). Recall, in addition, that 
the scattered-light intensity is directly proportional to the structure factor S(k) only in the case of unpolarized light 
scattering. Finally, recall (Sec. VI) that the values of the Porod amplitudes (but not those of the Porod exponents) 
in the general case of unequal elastic constants in the Frank free energy ( p6| ) are expected to differ from the values 
(derived in the one-constant approximation) listed in Tables || and |l| . 

For comparison, we also show in Tables | and || the results for the corresponding defects (when they exist) in 0(2) 
and 0(3) symmetric vector-model systems. In these systems, the normalization factor M norm is given by Vs 2 , where 
s is the order-parameter magnitude. Identical results for some of the O(N) defect configurations were previously 
obtained in Ref. Gal. 



TABLE I. Power-law contributions to the structure fa ctor arising from topological defects in systems with spatial dimension 
d = 2. Results are given per unit defect density (see Sec. VII). For biaxial nematics, Si and S2 denote the uniaxial and biaxial 
amplitudes of the nematic order parameter, and C x , C y and C z denote the three topologically distinct 180° defect types (see 
Sec. 0). 
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Biaxial 


0(2) vector 


0(3) vector 
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nematic 
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Uz ■ 8 3S 2 +S 2 K 


does not exist 


topologically 
unstable 


360° disclination point 


topologically 






topologically 




unstable 


fa ^ a \ ■ 8^2 (3S 2 -S 2 ) 2 , _ 4 
(Si < S 2 ) ■ 3 tt k 




unstable 
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TABLE II. Power-law contributions to the structure factor arising from topological defects i n sys tems with spatial dimension 
d — 3. The results given in the first three rows are given per unit defect density (see Sec. VII). The results in the fourth 
(fifth) row correspond to a single twist (wedge) circular disclination loop of radius 7? in a system of volume V; £ denotes the 
angle between the scattering wave-vector k and the loop axis. For biaxial nematics, 5*i and S2 denote the uniaxial and biaxial 



amplitudes of the nematic order parameter, and C x 
Sec. 0). 



C y and C z denote the three topologically distinct 180° defect types (see 
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We conclude this Paper by highlighting some features of our results. As discussed in Sec. [H], the value of the Porod 
exponent £ can be correctly anticipated purely on dimensional grounds. It should be noted, however, that the value 
of the dimensionless Porod amplitude A often differs considerably from unity. For example, in the case of the nematic 
hedgehog defect, we obtained A = 367T 4 ~ 3500. In addition, the Porod amplitude depends significantly on the order 
parameter in question (and not just on the spatial and defect dimensionalities). For example, the Porod amplitude 
for the ne matic hed gehog exceeds the amplitude for the 0(3) monopole in d — 3 by a factor of 37r. 

In Sees. [VE and [VF, we analyzed in detail the conditions under which the Porod tail of the structure factor is not 
overshadowed by the power-law contribution from transverse thermal fluctuations of the nematic director. [Similar 
considerations apply to any system possessing continuous symmetry of the order parameter.] We concluded that for 
experimentally accessible defect densities in a nematic phase-ordering experiment, the Porod tail should be observable 
over a range of 1 to 3 decades in k; the range of observability is limited by defect interactions at small k and by thermal 
fluctuations at large k. In the experiments reported in Refs. dll, ^ ne measured structure factor did exhibit a Porod 
tail over at least 1 decade in k, but the crossover at high k to the fluctuation-dominated regime was not observed. 
This was found to be consistent with our theoretical predictions under the conditions of these experiments. 

The crossover from S(k) cx k~ 5 (in a system where disclinations dominate over hedgehogs) to S(k) oc fc~ 2 (the 
thermal-fluctuation-dominated regime) was predicted in Sec. IV F to occur approximately at the scattering wave- 
vector value k u = (37r 3 /8A) 1 / 3 L^j 3 , where -Ldisci is the total disclination length per unit volume of the system. As 
this crossover is very sharp, and is not affected by defect interactions, measuring k u can give a precise estimate of 
the total disclination length. We are not aware of any use of this method in the experimental literature to date. 
The proportionality of the Porod tail to the defect density has been used previously |3]|| to extract the power law 
characterizing the decay of disclination length; extracting the absolute defect density in this way, however, would 
require the knowledge of the (rarely experimentally available) normalized structure factor. The location of the 
crossover k u , on the other hand, can be extracted from the unnormalized structure factor. 

In Sec. IVD, we contrasted the contributions to the structure factor arising from twist-type and wedge- type discli- 
nation loops. The wedge-type loop of radius R gives a k~ 6 contribution at k < R^ 1 in addition to the usual fc -5 
disclination contribution at k > R , as it has the structure of a nematic hedgehog at large distances from the 
loop. Likewise, any curved disclination segment of the wedge type gives rise to a fc~ 6 contribution. In Sec. IV F, we 
attempted to connect these findings to the observation that the structure factor in phase-ordering uniaxial nematics 
approaches the k~ 5 power law at high k through effective exponents larger than 5 (and close to 6) at intermediate 
k. Such a behavior of the structure factor had been found to occur in experiments as well as in numerical 
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simulations |27|] , but not in approximate analytical theories M. We concluded that it was, in principle, possible, but 
nevertheless unlikely, that the "approach from above" to the exponent x = 5 observed in the experiments was due to 
k~ e contributions from hedgehog defects. It remains a challenge for future work to determine whether this behavior 
of the structure factor is caused by hedgehog defects, curved wedge-type disclination segments, or correlations that 
fall beyond the range of the Porod regime. 
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